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1 Introduction 



The standard sandpile model introduced in |l| became of mathematical interest 
after the paper by Dhar who discovered its Abelian structure. The model is 
defined on a finite hypercubic ci-dimensional lattice. Each site i is characterized by 
a nonnegative integer variable Zi called the "height". If Zi < 2d for all i, the sandpile 
is said to be stable. A vertex is picked at random and its height is increased by one. 
If Zi > 2d, then the site is unstable and topples giving one particle to each of 
its neighbors which in their turn, can be unstable and topple. The process called 
" avalanche" continues until all sites become stable and then a new particle is added 
to the lattice. If the avalanche reaches the boundary, sand disappears in a sink 
connected with all boundary sites. 

The "size" of an avalanche may be measured by the total number of topplings 
s, the number of distinct sites toppled a, the diameter of the region affected by 
avalanche r, and duration of avalanche t. It is generally believed, the probability that 
the avalanche has size x{x = s, a, r, t) varies asymptotically as x~^(r = r^, r^, r^, r^) 
when x — >■ cx). The mean- field value of the exponent obtained from exact solutions 
on the Bethe lattice 0] and on the full graph |^ is 3/2. It is expected that = 3/2 
also for Z'^ when d > d^ where d^ is the upper critical dimension. 

The first attempt to find du was made by Obukhov soon after the sandpile 
model was proposed. Briefly, Obukhov's arguments can be summarized as follows. 

(i) Consider changes resulting from an avalanche propagating through the lattice. 
If the system is in the recurrent state, the average sensitivity to a new excitation does 
not change. During the avalanche, the sites which have already toppled have heights 
lower than average, whereas the heights in the neighboring sites are larger than 
average. Hence, the previously activated sites repulse the new activation process. 
This situation resembles the True Self- A voiding Walk (TSAW). 

(ii) The essential difference between the sandpile and TSAW is a possibility of 
branching the activation process. Thus, the avalanche can be pictured as a self- 
avoiding branching process. 

Using the renormalization group and e-expansion, Obukhov calculated one- 
loop corrections to the mean-field theory and came to the value du = 4. This 
conclusion was supported by Dias-Guilera who analyzed non-linear stochastic 
differential equations derived from the models with continuously distributed heights 
[0]. Later on, Christensen and Olami [g] suggested du = 6 from an analogy be- 
tween spreading of avalanches and percolation. The mean-field treatment of a self- 
organized branching process was discussed in 0. 

The proof of correspondence between the branching TSAW and an avalanche 
process in the sandpile model meets two considerable difficulties. First, a typi- 
cal avalanche contains multiple topplings which violate self-avoidance of branches. 
In higher dimensions, multiple topplings are suppressed |0,|jl5|. However, then, 
the second problem arises. The branching process corresponding to an avalanche 
depends deterministically on a recurrent configuration where it was initiated. All 
recurrent states in the Abelian sandpile model have equal probabilities 0]. If self- 
avoiding branching processes having equal number of steps, have equal statistical 
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weights, they can be considered as lattice trees also known as branched polymers. 
The upper critical dimension of branched polymers = 8 was predicted in and 
then rigorously determined in pJ| , [p!2| , p!3[ . 



Thus, avalanches either do not correspond to the branching TSAW due to mul- 
tiple topplings, or do not have the specific TSAW weights due to equal probabilities 
of recurrent states. 

The computation situation is not less controversial. For the most part of numeri- 
cal experiments, the accuracy is not sufficient to distinguish between non-trivial and 
mean-field values of critical exponents in high dimensions. Only high-statistics data 
on large lattices provide some information. Grassberger and Manna [|14| investigated 
critical exponents for li < 5 and concluded that du = 4. Liibeck and Usadel [|l^ 
found du = 4, mentioning that the largest considered size of the system {L = 80) 
for d = 4 is too small and the corresponding avalanche distribution exhibits a very 
narrow power-law interval. In his paper |TB[, Liibeck pointed an important role of 
logarithmic corrections to the scaling at du = 4. 

At the same time, Chessa et al. performed a numerical study of critical 
exponents in dimensionality ranging from d = 2 to d = 6 and observed the mean- 
field behavior only in ci = 6 excluding, therefore, d^ = 4:. 

Very recently, Vespignani et al [jT8[ have derived du = 4 from a phenomenological 
field theory, refiecting the symmetries and conservation laws of sandpiles. 

An apparent inconsistency between the simple self-avoiding branching process 
and real avalanches, as well as contradictory numerical results of different groups, 
set to find a more transparent proof of the upper critical dimension for sandpiles. 

In this paper, we prove the upper critical dimension d^ = 4 using Lawler's the- 
orems for intersection probabilities of random walks and loop-erased random 
walks. We introduce again the self-avoiding branched polymers for description of 
avalanches and show that avalanches are spanning subtrees embedded into a span- 
ning tree of the whole lattice rather than usual lattice trees. The problem of fractal 
dimension of avalanches is reformulated as that for the spanning subtrees. Using 
Majumdar's result of the equivalence between the chemical path on a spanning 
tree and the loop-erased random walk [^, we reduce the problem to estimations of 
intersection probability between random walks and loop-erased random walks. 



2 The Model 

We consider the sandpile model on the d-dimensional hypercube A C Z'^. Elements 
of the state space {0, 1,2, ...{2d — 1)}^ are called stable configurations and are de- 
noted by C. The value C{i) = G A is the height of the sandpile at the site i. 
Given a configuration C and a lattice site i, aiC is the stable configuration obtained 
by adding a particle at i, and relaxing the system by topplings at all unstable sites 
j, Zj > 2d. On toppling at the site j, 

Zi^ Zi- Aij i e A (1) 
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Elements of the matrix A are: An = 2d for all i G A; Aij = —1 for all bonds 

The operators commute with each other 

[ai,aj] = i,jeA (2) 
This allows one to define the identity operator 

J, = n «f " ^ e A (3) 

j€A 

and the equivalence relation between two configurations C' and C" 

C'{j) = C"{j) + J2niA,, 3 eh (4) 

i 

where G A are integers. The equivalence means that both C' and C" tend to 
the same stable configuration after topplings at all unstable sites. 

The set of stable configurations {C} splits into two subsets: {C} = {C}R[j{C}T 
where {C}ji is the recurrent set and {C}t is the transient set of the sandpile process. 
In {C}r, one can define an inverse operator a^^,i G A. Then, the operators {ai} 
form a finite Abelian group [Q]. 

Due to Eq.(^), each C G {C}r is an element of the super-lattice in the |A|- 
dimensional Euclidean space. The basis-vectors of this lattice are the rows of the 
matrix A. Therefore, the number of elements in {C}r is 

Nr = DetA (5) 

The recurrent configurations in {C}r have equal probabilities N^^. 

Denote by dA the set of boundary sites of A and by B the set of elementary 
cubes of the (i-dimensional Euclidian space centered at i G dA. We choose a site 
i G Z*^ not belonging to A and call it the root. The faces of cubes from B which 
can be connected with i without intersections with another face, are called external 
faces. We connect each site i G dA with i by z/j bonds where z/j is the number of 
external faces of the cube i. The lattice A together with i and new bonds forms a 
graph denoted by G. The site i is a sink for particles leaving the lattice A during 
the avalanche. 

Using the explicit form of the identity operator Eq.(|^), one can construct a new 
identity operator 

= n (6) 

IdA says that adding Vi particles to each boundary site triggers an avalanche of a 
special form: each site i G A topples exactly once and the initial configuration C 
remains unchanged. One can use Iq\ to construct a graph representation for any C G 
{C}r by the so-called "burning algorithm" |]22|. Consider the topplings initiated by 
JgA as a fire starting at i and burning sequentially all lattice sites. Once the rules 
of propagation of fire are fixed, the set of bonds along which fire propagates forms a 
spanning tree T of the graph G. There is one-to-one correspondence between {C}r 
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and the set of spanning trees T. Then, Dhar's formula, Eq.(^, for coincides 
with Kirhhoff 's theorem for the number of spanning trees . 

The matrix A in Eq. is the minor of the Laplacian matrix of the graph G 
corresponding to the element (Ac)^^. According to Dhar 0, the Green function 

G^, = [A-% (7) 

is the average number of topplings at the site j due to a single particle added at i in 
a configuration C G {C}r. For the random walk defined by the matrix A, Gij is the 
expected number of times the walk started at i visits the site j until it is trapped 



by the absorbing site i [26 



To find the spanning tree representation for G^j, we delete from A the elementary 
cube centered at i and consider its 2d faces as a part of the new external boundary. 
The sites adjacent to i form a new boundary set d^A. Repeating the construction of 
a spanning tree by the identity operator, Eq.(^, with the boundary set dA U d^A, 
we obtain a two-component spanning tree having the roots i and i for two different 
subtrees. 

Proposition 1 proved in (see also |2^) reads: For any connected graph G 
with a fixed vertex i 

G,, = Ar(^'^-)/l^l (8) 

where A^^*'-?) is the number of two-component spanning trees having the roots i and 
2, such that both vertices i and j belong to the same component; |T| is the number 
of spanning trees on G. 



3 Avalanches and Waves 

An avalanche starting at i is the process of transformation C aiC. Generally, 
an avalanche consists of multiple topplings at different sites and has a complicated 
structure. However, one can try to decompose it into simpler subprocesses. Suppose 
Zi = 2d after adding one particle to i. Since the topplings can be performed in any 
order, we topple once at i, and then topple all other unstable sites keeping the site 
i out of the next toppling. This sub-avalanche is the first wave of topplings. If i 
is still unstable, we topple i once again giving rise to the second wave of topplings. 
This process is continued until i becomes stable. Thus, an avalanche is broken into 
a sequence of waves of topplings. 

Three important properties of waves make them useful for the analysis of avalanche 
statistics [^. Let 5*^, > 1 be a set of sites toppled during the k-th wave. Then, 

(i) each j G Sk topples once and only once; 

(ii) both 5"*; and A \ Sk are connected sets; 

(iii) k-th wave is the last wave of an avalanche iff the initial site i & Sk has a 
neighbor in A \ S'fc. 

The number of waves in the avalanche C — > OjC will be denoted by ni{C) or 
simply rii. 
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The spanning tree representation of waves can be constructed as follows. Given 
a configuration C G {C}/j, PiiC) represents the projection of C on the lattice A \ i 
where the site i is now considered as the second root. Let Ck, k < riihe an (unstable) 
configuration obtained after the k-th wave starting at i is completed. Then, Pi{Ck) 
can be produced from Pi{C) by the operator 



W, = {U a,)' (9) 



acting on the lattice A\i. 
Consider the configuration 



P,(Cfc_i) = Wk-lP^{C) (10) 

preceeding the k-th wave. If one applies to Pi(Cfe_i) the identity operator 

IdAud^A = ( n n (11) 
jed^A jedA 

one obtains a two-component spanning tree representing Pj(Cjt_i). On the other 
hand, the expression in the first brackets is exactly the operator providing the k-th 
wave. Therefore, the subtree rooted at i and embedded into the subtree rooted at 
i, is the graph representation of the k-th wave. 

This result gives another proof of Eq.(|]). As each wave beginning at i and 
involving j corresponds to one toppling at j, N^'^'^y\T\ is the average number of 
topplings at the site j for a single particle added at i. Then, Eq.(||) follows from 
Eq.(0). 

Due to the graph representation, we come again to a lattice-tree or a branched- 
polymer picture of avalanches. However, the essential difference with the qualitative 
Obuhkov arguments is, firstly, that we associate the lattice trees with waves but not 
with the whole avalanches, and secondly, the lattice trees are conditioned by the 
spanning tree construction. 

Let Sk denote the number of sites in the k-th wave Sk- The total number of 
topplings s in an avalanche is 

rii 

s = Y.^k (12) 

k=l 

where n is the number of topplings at the initial point i or, equivalently, the num- 
ber of waves in the avalanche. It is expected that, at least for sufficiently large 
dimensions d, the probability distribution of s has the scaling form 

Pis;L)r~.s'^fis/L^) (13) 



where D is the capacity fractal dimension of avalanches. It follows from Eq.(|T2|) 
that D depends on the fractal dimension of waves and the number of waves in an 
avalanche. In c? = 2, the dimension of waves is d,,, = 2 and the number of waves is 



expected to obey the scaling law ||22| 



(ri.) ~ (14) 
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Then the exponents Ts,Ta,Tr,Tt can be determined in terms of y. In higher dimen- 
sions, multiple toppling events occurs more rarely showing a tendency to decrease y. 
In Section 7, we will show that for d — i, {ui) grows with r not faster than logarith- 
mically. If so, the problem of critical dimension for avalanches can be formulated as 
that for waves. 



4 Statistics of Waves 

We start this section with the definition of the fractal dimension of waves. Consider 
the graph representation of a wave as a two-component spanning tree on the graph 
G with the roots i E A and i. The wave S is the set of s = l^l sites belonging to 
the component rooted at i. The radius of the wave is 

R{S)^sup{d{i,j):j eS} (15) 

where d{i,j) is the distance between the sites i and j. We will say that waves have 
a fractal growth if 

(s) ~ (16) 

for large R, where dfg is the dimension of fractal growth and the average is taken 
over all waves of the radius R. The fracral growth, however, does not guarantee the 
fractal structure by itself. 

Surround the point i by a ball of radius r and consider all sites of S inside the 
ball. Denote the number of internal points by m(r). We say that waves have the 
fractal density if, given e > 0, there exist r(e), R{e) and a constant c < 1 such that 

\log{m{r)) - dfjogr] < e (17) 

on the interval r(e) < r < cR for all R > R{e). 

If dfg = dfd = df, the fractal is uniform and df is the fractal dimension of 
waves on the (i-dimensional lattice. We do not prove here uniformity of the fractal 
structure of waves. Assuming it, we will find the upper and lower bounds for df. 

In the spherically symmetrical case, we can introduce the density of waves which 
varies as 

in the scaling interval r(e) < r < cR. 

If df is known, one can determine the probability distribution of waves P{x),x = 
s,a,R,t {s = a for waves). Indeed, the expected number of topplings at a point j 
in an avalanche starting at i, |i — j| = r varies as 

(n(r)) -p(r) H P{x)dx (19) 

Jr/c 
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where the integral gives all waves having reached the radius r/c. On the other hand, 
(n(r)) coincides with the Green function, Eq.(|^). The asymptotics of G{r) in the 
(i-dimensional case is 

Gir) - ^ (20) 



Comparing with Eg . ([191) , we find that 



PiR) ~ (21) 

The probability distribution for s follows from Eq.(|l6D provided that s ~ when 
R —>■ oo 

Pis) -L (22) 

with Tg = 2 — 2/df. The exponent reaches its mean-field value 3/2 when df = 4. 

To proceed with the determination of df, we need a more elaborate definition 
of p(r). Consider a wave S with the initial site i. The set of sites not belonging 
to S is denoted by S, S U S = A. According to Eg. (]T3|) , a wave has the radius 
not less than R if there exists a site i & S , \i — i \ = R. The wave S has the 
density less than 1 at the radius r < i? if there exists a site i E S, \i — i \ = r. 
Consider the two-component spanning tree with the components corresponding to 
S and S. The path T{i,j) on a tree between the points i and j is the seguence of 
bonds {ii,i2), ...{in, j)- Any two points belonging to one component can be 

connected by a unigue path. Then, we can define the density of waves at the radius 
r < R hy the conditional probability that, given a path T(i,i ), |i — i | = R, there 
exists a path T{i,i ), \i — i \ = r, 

p{r) = 1 - Prob{T{iJ)\T{i,i')) (23) 



The expression in the right-hand side of Eg. ([23| ) is the non-intersection probability 
of two self-avoiding paths on a tree, called often "chemical paths". The first path 
T{i,i ) connects the site i with the site i belonging to the wave S; the second path 
T{i,i ) connects the site i not belonging to 5* with the sink i. Having no rigorous 
results for this probability, we can, however, reduce it to non-intersection probability 
between one chemical path and the simple random walk. 

For this purpose, we need a generalization of Proposition 1 to the multicom- 
ponent case. Consider a connected graph G and select a subset of its vertices 
{v} = ii,i2, ...,iu. Let ^{v} be the matrix obtained from the Laplacian matrix 
by deleting all rows and columns corresponding to ii,i2, ...,iv. The Green function 

= l^^h (24) 

is the expected number of visits of the site j by the random walk starting at i before 
it is absorbed at one of the sites ii,i2, ...,iu. 

Proposition 2 For any connected graph G with a fixed subset of vertices {v} 

G.J = (25) 
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where N^J^ is the number of {u + l)-component spanning trees having the roots 
ii,22, ■■■,iu and i such that vertices i and j belong to the same component; |T{t,}| is 
the number of i/-component spanning trees on G having the roots ii, 12, ■■■,iu- 

Remark If a site G {f } is an isolated site of the spanning tree, we consider it 
as a component consisting of the single site. 

Proof Consider the graph G which is the graph G with all sites ii,i2,...,iu 
contracted to the single site iq. Proposition 1 is applicable to G . It relates the 
Green function, Eq.(p4D, with the number of two-component spanning trees having 
the roots io and i, such that the sites i and j belong to the same component. 
Returning to the graph G, we get from each configuration on G , a spanning tree 
where the component containing iq splits into u components. This construction 
implies Proposition 2. 

In order to treat the right-hand side of Eq. (|23| ) as a random walk probability, we 
consider the graph G with the subset of vertices {v} coinciding with the set of all 
sites of the chemical path r{n ). Then, by Proposition 2, we have 

G,i = (26) 

where G-^'i is the expected number of visits of the site i by the random walk starting 
at i and escaping the chemical path r(u ). If z/ is the number of sites in r{ii ), then 

A^l*l is the number of + 1) component spanning trees having sites i and i in one 
component. The remaining u components can be joint into a single component if 
one adds to u spanning subtrees rooted at the sites of T{n ), the path T{ii ) itself. 

Then, N^^y' is the number of two-component spanning trees which have the path 

r(ii') in one component corresponding to the wave S and the path T{n) in the 
second component S having the root i or, equivalently, i. \T^v}\ is the number of 
all one-component spanning trees containing the path r{ii ). 

Along with |T{„}|, we introduce N^^}, the number of two-component spanning 

trees containing the path r(u') in one component and the root i in the second one. 
By definition, 

Prob{T(vi'm^^')) = ^ = (27) 

Prob{T(ii')\Tiit)) = ^ (28) 

where the Green function G^ is the expected number of returns to the initial point 
for random walks starting at i and escaping the path r{ii ). 



Using Eq.(^, we get 



Finally, we note that the ratio of Green functions in Eq.(28) is the probability 
that the random walk starting at i escapes the path r{ii ) and reaches the site i. 
This fact follows from the property of Green functions known for infinite lattices 
as the "weak ergodic theorem" |2^. Denoting hj F{i the escaping probability 



averaged over all paths r(z,i'), we can summarize the obtained results as 
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Proposition 3 The waves starting at a site i of the lattice A with the sink i 
have the density 

p{r) = 1- F(i'\ii') (29) 
at the point i' , \i — i'\ = r, if the radius of the wave exceeds R = \i — i'\ > r. 



5 Loop erased random walks 



The loop erased random walk (LERW) has been introduced by Lawler [21] to modify 
the uniform measure of the usual self-avoiding random walk where every possible 
walk of a fixed length is given the same statistical weight. Let [xQ,...,Xm] be the 
simple random walk in Z'^. Then the LERW can be constructed as follows: let j be 
the smallest value such that xj = Xi for < i < j < m. Deleting all steps between i 
and j, we obtain a new walk [xq, ...,Xi, Xj+i, Xm]- If the new walk is self-avoiding 
we stop the process; otherwise, we continue the loop-erasing operation until we get 
a self-avoiding path. 

The distribution for each step of the LERW depends on the entire past history 
being equivalent to that given by the transition probability of the Laplacian random 
walk |27|,[^]. Another equivalence has been found by Majumdar |2^ who consid- 



ered a model of growing trees introduced by Broder . The chemical path on these 
trees corresponds to the LERW. Like the burning algorithm, Broder's algorithm in 
the limit of lattice filling gives spanning trees with equal probability. Therefore, the 
statistical properties of the LERW are identical to those of the chemical path on 
spanning trees. 

Let Rn{i) and Ln{i) be simple and loop-erased random walks starting at the site 
i and let Hj and Fj denote paths of the walks 

Ui{a) = {Rn(i) : < n < a} (30) 

Ti{a) = {Ln{^ : < n < a} (31) 

Denote by (y{n) the number of steps of the simple random walk needed to produce 
n steps of the LERW 

a{n) = sup{j : L„(ii) = Rj{i)} (32) 

For any random walk, (y{n) > n. For a fixed cr(n), Ili{a{n)) is the path of a finite 
random walk. Then, Fj(n) is the finite LERW obtained from Ili{cr{n)) by the loop- 
erasing procedure. 

The escaping probability of a set A G Z'^ is defined as 

Es{i, A) = Prob{Rr,{i) ^ A, n = 1, 2, 3, ...} (33) 

where i is the starting point of the random walk. Putting A = Fj(n), we can define 

Q^j^i as 

an = E[Esit,T,{n))] (34) 
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i.e. a„ is the probability that an infinite independent simple random walk starting 
in z G Z"^ does not intersect the finite LERW derived by erasing loops from cr(n) 
steps of an independent random walk also starting in i. 

For d = 4, Lawler |^ has found the lower and upper bounds for a„ 



lim mf- — ; > — (35) 

n^oo log log n 2 

and 

r log . 1 

hm sup- — < — - [Sbj 

n^oo log log n 3 



The upper bound Eg . (|36[) was conjectured to coincide with the exact asymp- 



totics. This result would imply that the mean square distance for LERW's behaves 
as 

(0~n(logn)V3 (37) 

whereas (r^) ~ n for d > 5. 

It is convenient to bring here the known intersection probability for two simple 



random walks ||T9|. Let Rn{i) and Rnij) be independent random walks starting at 
the points i and j, \i — j\ = r. Then, there exist constants Ci > and C2 > such 
that 

ci < Proh{T[i{n) n Ilj{n) ^ 0} < Cs (38) 

for ci < 4 

ci(logn)"^ < Proh{Ili{n) n Ilj{n) ^ 0} < C2(logn)"^ (39) 

for ci = 4 

^^^(4-d)/2 < Proh{Iii{n) n n,(n) ^ 0} < csn^^-'^)/^ (4Q) 

for ci > 4, if a^/n < \r\ < by/n, where a and b are positive constants. 

The non-intersection probability of a finite random walk and an infinite random 



walk starting at the same point i for ci = 4 is 1 19 1,| 32 



Prob{Ui{n) n ni(oo) = 0} ~ c(log n)"^/^ (41) 

Also, for further estimations, we need the intersection probabilities between a finite 
random walk starting at i and an infinite random walk starting at j for ci = 4 |[3l[ 



Prob{Ui{n)nUj{oo) 7^0} < Ci log(l + l/ct)/ logn (42) 
if a = \i — j\'^/n. 



6 Upper and lower bounds for density of waves 

Every random walk on a graph G with an absorbing set {v} and the sink i is trapped 
either by {v} or by i. Then, Proposition 3 means that p(r) is the intersection 
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probability between the random walk starting at the site i ,\i — i\ = r and the 
chemical path r{i, i ) , \i — i \ = R. 

Majumdar I^O] has proved that chemical paths correspond to LERW's which, 
in turn, can be obtained from simple random walks by the loop-erasing procedure. 
Let Rn{i) be the random walk corresponding to the LERW Ln{i)- As the number 
of steps a{n) in the path Ili{a{n)), Eg . (|30|) , needed to reach the radius R, exceeds 
the number of steps n in the LERW path Ti{n), Eg. (PT]) , we can estimate p(r) as 

p(r) < Prob{Ui{a{n)) n (oo) ^ 0} (43) 

In d = 4, we get from Eg. (42) 

p(r) < c log(l + -)/ log r (44) 

a 

where a = r'^/a{n) ~ r'^/R^. The fractal density of waves in d = 4 decays with 
r at least logarithmically when a is fixed. In > 4, we can use Eg.(^OD assuming 
that the asymptotics does not change when one of the random walks is extended to 
infinity. Then, we have 

p(r) < cr^-^ (45) 

From the definition, Eg.(|TB|), we can see that the fractal dimension of waves df < 4 
for all (i > 4. 

In d < A, the intersection probability. Eg. (|55D , and, therefore p(r), are restricted 
from above by a constant. 

To get a lower bound for p(r) in = 4, we will use the upper bound Eg. (|36|) 
for escaping probability a^. Consider the random walk Rn{i) and decompose its 
path nj(oo) into two parts: nj(oo) = Ili{k) U nj'(oo) where k is the moment of the 
first hitting into the site i separated from i by the distance r: Ri{k) = i ,Ri{l) ^ 
i ,1 = 1,2, k — 1. Then, escaping probability a„ for the LERW of length n by the 
random walk Rn{i) is the product of escaping probabilities an{i,i ) and a„(i , oo) 
before and after the first hitting into the site i 

a„ = a„(i,/)a„(z', oo) (46) 

Granting that a„(i', oo) = F(i'\ii'), we can write Eg. (p9D as 

P{r) = 1 - (47) 

where n is the length of the LERW coinciding with the path T(i,i'). 

The probability a„(2, i ) to reach the point i from i for m steps, m ~ I? — i p = 
r^, escaping the LERW of length n is not less than the probability to avoid the 
random walk of length cr(n) from which the LERW was obtained by the loop- erasing 
procedure, and therefore not less than the probability Prob{Ili{oo) fl Ili{m) = 0} to 
avoid an infinite random walk. Using Egs. (|36|) and (^T]) we get from Eg. (|47|) 

P(r) > 1 - e 

(log 77.)^/'^ 
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For any r ~ m^^^, we see from Eq.(|37|) that the number of steps n in the LERW 
increases with the radius of wave as (with the logarithmic correction) and p{r) 
approaches 1 when R ^ oo. The only fractal dimension which is consistent with 
both upper and lower bounds Eq.(^^ and Eq.(^8D is 4 although the fractal density 
Eq.(|T8|) as well as the probability distribution of waves Eqs.(^),(^) need logarith- 
mic corrections. 

The lower bound Eq.(^) becomes stronger in d < 4 because of decreasing es- 
caping probability F{i \ii ). Assuming that dj is a non-decreasing function of d, we 
conclude that the fractal dimension of waves is 

df = 1^ for d < 4 ^ ^ 

^ 14 fore/ > 4 ^ ^ 

According to Eq. (p2D , the upper critical dimension for waves is du = 4. 



7 Number of Waves in Avalanches 

In order to derive the upper critical dimension for avalanches from that for waves, 
we need some information about the expected number of waves in an avalanche. 
Due to the property (iii) of waves mentioned in Section 3, the probability that a 
wave starting at the site i is the last wave in an avalanche is proportional to the 
density of "holes" in the wave at its origin, i.e. the probability that a site in the 
vicinity of i does not belong to the wave. By Eq.(^), this probability is F{i \ii ) 
which coincides simply with the escaping probability a„ in the case r = \i — i \ ^0. 

Avalanches of radius r consist of the waves whose radius does not exceed r. 
Therefore, the expected number of waves can be estimated as 

{rii) < ca-^ (50) 

where ~ ?T,(logn)^/^. Using Eq. (|55D , we get 

{rii) < c(logn)^/^ (51) 

or 

W <Ci(logr)i/2 (52) 

Thus, in d = 4, the number of waves in an avalanche of radius r grows not faster 
than logarithmically. 



In d = 2, the distribution of last waves is known |33 



^K^) ~ ^ (53) 

The probability for the wave of radius R to be the last in an avalanche is Pi{R) 
divided by the general distribution of waves Eq.(^). Then, the expected number of 
waves (rii) in d = 2 has the upper estimate 

{m) < cr'^l^ (54) 
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which is consistent with Eq.(14) where the exact value y = 1/2 was conjectured |34 

For (i > 4, the escaping probabihty a„ remains finite for all n. Therefore, the 
expected number of waves in an avalanche is restricted by a constant. 

The logarithmic growth of (nj) in c? = 4 given by Eq.(|52D implies that the proba- 
bility distributions for waves and avalanches can differ not more than by logarithmic 
correction. 

To estimate this correction, we take instead of the density p{r) its upper bound 
Eq.(|4^) if r.f.s. Eq.(^4[)is less than 1, and put p(r) = 1 otherwise. Comparing Eq.(|20D 
with Eq.(|l9]) where the new p(r) is taken, we get the probability distribution for 
waves in the form 

PW = (55) 

Also, we take instead of the lower bound Eq.(p5D for a„ its exact asymptotics con- 
jectured by Lawler |30|. Then, 

{n,) ~ (logr)^/3 

The maximal difference between avalanche and wave distributions corresponds to a 
situation when all waves in an avalanche of radius r have the maximal radius r 

{ni)Pavai{r)dr = P{r)dr (57) 

Hence, the avalanche distribution in d = 4 can be estimated as 

PaUr) ~ (58) 

with 2/3 < 7 < 1. 

The probability distribution of the total number of topplings in a wave P{s) 
follows from Eq.([55|) provided the leading asymptotic of s is s ~ R^/ logR: 

(logs)^/2 

P(s) - (59) 

In > 4, the asymptotics of the avalanche and wave distributions coincide and 
correspond to the mean- field behavior with the exponent Ts = 3/2. 
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